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Abstract 

We study the three-dimensional quantum Ising spin glass in a transverse 
magnetic field following the evolution of the bond probability distribution 
under Renormalisation Group transformations. The phase diagram (critical 
temperature Tc vs transverse field V) we obtain shows a finite slope near T = 0, 
in contrast with the infinite slope for the pure case. Our results compare very 
well with the experimental data recently obtained for the dipolar Ising spin 
glass LiHoo.i67Yo.833F45 in a transverse field. This indicates that this system 
is more apropriately described by a model with short range interactions than 
by an equivalent Sherrington-Kirkpatrick model in a transverse field. 

PACS: 75.50.Lk, 05.30.-d, 75.10.Jm, 75.30.Kz 
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The role of quantum fluctuations in spin glasses has been a long standing theoreti- 
cal problemEhi. The so-called proton glasse^^ — a random mixture of ferroelectric and 
antiferroelectric materials such as Rbi_2;(NH4)2;H2P04 — provide an experimental realiza- 
tion for quantum spin glasses. Within a pseudospin description of such hydrogen-bonded 
systems, the proton position in the two potential minima is represented by Ising states, 

= ±1, and the tunneling between the minima by a transverse fleld term, Ta^, where F 
is the tunneling frequency^. The theoretical study of Ising spin glasses in a transverse fleld 
(TISG) has then attracted renewed interest, especially in relation to the analogue of the 
Sherrington-Kirkpatrick modeli in a transverse fleld (TSK); see, e.g., Ref. ^ and references 
therein. More recently, the magnetic susceptibility of the (long-ranged dipolar) Ising spin 
glass LiIIoo.i67Yo.833F4 has been measured in the_presence of a transverse fleld Ht, from 
which a phase diagram Tc{Ht) was determinedSo. Therefore, it is of interest to discuss 
the main differences between the phase diagrams of the transverse Ising model in both pure 
and spin-glass cases. In view of the long range nature of the interactions in the dipolar 
glass, we are also particularly interested in establishing whether this system can be suitably 
described by a short-range model or one has to resort to the TSK model. Here we address 
these questions using real-space scaling methods. 

The TISG model is described by the Hamiltonian: 

TV N 

-m = Y.J.,ata^^+VY.at (1) 

i<j i 

where the , n = x, z are Pauli spin matrices, F is the transverse fleld, i and j are nearest- 
neighbor sites on a simple cubic lattice, and the Jij are uncorrelated exchange couplings 
chosen at random from an even distribution. For zero transverse fleld the model reduces to 
the classical Ising spin glass, and quantum effects are brought in by increasing the fleld. At 
flnite temperatures, the effect of the transverse fleld is to depress the spin glass transition 
temperature, whereas at zero temperature quantum fluctuations are the only mechanism 
driving the system to a phase transition at a critical value of the transverse fleld. It is 
interesting to note that the lower critical dimension (di) for zero temperature transitions in 
the Transverse Ising spin-glass model igil c?^ = 1, unlike the classical Ising case 0, df = 3. 

In the context of real-space renormalisation group, the simple cubic lattice may be ap- 
proximated by hierarchical Migdal-Kadanoff ce llsEl; see Fig. 1. The terminal sites are con- 
nected by b'^~^ bonds 'in parallel', each of which consists of b bonds 'in series'; b is the scaling 
factor (6 = 2 in Fig. 1) and d is the space dimensionality of the lattice (c? = 3 in this case). 
For the transverse Ising model (TIM), the non-commutation aspects are present at the clus- 
ter level, in the sense that individual spins are not in a deflnite state. This can be dealt 
with by referring the density matrix to the basis |mim2 . . .mN), where cr||mj) = mj|mj), 
and deflning the renormalisation group transformation (RGT) by the mapping of diagonal 
elements only. This approach has been successfully used in a detailed study of the pure 
and bond-diluted TIM in two dimensionsi^. In the present work bond disorder is included 
within a statistical renormalization group (SRG) treatment: one follows the effect of a RGT 
on the probability distributions of the relevant parameters, instead of forcing them into a 
particular formfij. Several aspects of the Ising spin glassH and of the random fleld Ising 
mo del0 have been elucidated by treating disorder this way. 
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Thus, for a given bond configuration {{Jij}), an RGT for the d = 3 system (Fig. 1) is 
defined by 

(mim6|p'(K')|mime) = (mime|p(K)|mim6) (2) 

where K' = {J',r',C') are the renormahzed quantities in the two-site cell, K = ({Jjj},r) 
refers to the original cluster, and 

(mim6|p(K)|mim6) = ^ {mim2msm4m^mQ\p{Ji)\mim2rnsrn4'mr,rnQ) (3) 

m2m3m4m5 

is obtained by performing the partial trace on the internal spins, keeping those on the 
terminal sites fixed (see Fig. 1). The third matching condition required to solve the problem 
for the three unknows (J', F', C) is obtained by preserving the thermodynamical average of 
the rotation operator 

For the renormalized cell, Eqs. (^ and (^) provide analytical expressions for the three 
primed unknown variables^. Since disorder destroys the point group symmetry of the orig- 
inal cluster, the Hamiltonian matrix is written in a 64 x 64 representation and Eqs. (|^) and 
(I) are calculated numerically. At this point we should comment that a single renormalized 
field appears in the recursion relations as a result of the approximation employed here: the 
fields are assumed to be uniform in both renormalized and original cells. This assumption 
can be justified, to some extent, by recalling that the transverse magnetic field behaves as 
an irrelevant variable in the pure TIM. Since the transverse field is not a symmetry break- 
ing operator, this should hold in the spin glass case as well. On the other hand, we could 
have allowed the field at each site to follow an RG trajectory. The irrelevance of the field 
would manifest itself through a distribution evolving from an initial delta function centered 
at F = Fq into one centered at F = 0. 

The initial probability distribution for the exchange couplings is 



P{Ji) = . ^ exp 

V2^ 



J? 



2J2 



(5) 



where we have assigned an index i to each bond in the original cell and taken ( Jj J;) = J'^Su 
and (Jj) = 0. We start the iteration by choosing eight bonds distributed according to 
Eq. and one value of the transverse field (-P(F) = 6(Ti — F)), to feed the recursion 
relations and generate a new value of the field and of the exchange coupling. This procedure 
is repeated about 10,000 times and we obtain two renormalized distributions, -P'(Jj') and 
P'(F'). We make use of these distributions to feed the recursion relations in the next step 
of the renormalization process: eight bonds and one field are chosen according to the new 
distributions, and this is done again 10,000 times. The evolution of the distributions is then 
followed along the renormalization process. We find that the bond probability distribution 
remains symmetric around J = at each iteration; that is, (Jij) = 0. The attractors of the 
different phases are determined by fixed distributions characterized by their width J and 
mean value (F), as follows: 
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J — >■ oo, (r) — > spin glass 

J — i> 0, (r) paramagnetic. 

so that the critical curve is determined as the boundary between these two different be- 
haviours. Note that in every case the average value of the transverse field distribution 
iterates to zero. 

In Fig. 2 we present our results for the phase diagram, vs T, for both pure and spin 
glass cases. For comparison we also display the spin glass data for the replica-symmetry- 
breaking (RSB) solution of the TSK modeli, and for LiHoo.i67Yo.833F4li^l0- In analysing the 
experimental data one should have in mind0'0 that the applied transverse field Hf gives 
rise to a level splitting ~ Hf (at low fields) which, in the context of Eq. (|l]), is proportional 
to F; thus, F ~ Hf. For the pure TIM, the critical line both for small fields and near 
T = has a square fit: T'c(O) — T'c(F) ~ F^, and Fc(T) — Fc(0) ~ T^, respectively. Still for 
the pure case, the critical field for zero temperature transitions obtained with the present 
RG is (F/J)c = 3.40, which should be compared with the series resultjil, (F/J)^ = 5.14; 
similarly, our result Tc(0) = 3.83J should be compared with Tc = 1.13J, obtained from 
series expansions for the three-dimensional Ising mo deli. As usual, the critical parameters 
obtained within a simple Migdal-Kadanoff approximation are quite inaccurate, but one is 
generally able to describe the qualitative aspects of phase diagramJlUli. 

For the transverse Ising spin glass, we obtain a curve Tc(F) with a finite slope near T = 0, 
unlike both the pure case and the replica symmetry-breaking-solution to the TSK model; 
see Fig. 2. Overall, the experimental data are better represented by the present approach 
than by the infinite-range mean-field mo deli. This can be explained by the fact that in the 
actual crystal, the dipolar interactions fall off with the distance, being effectively reduced to 
zero for distances greater than a few lattice spacings. In contrast, the interactions between 
any pair of spins in the mean-field model have the same intensity, irrespective of the distance 
between them. The calculated critical parameters in this case are: Fc = 1.58 J at T = 0. 
and Tc = 0.884J at zero transverse field. In agreement with the experimental resultsiyO 
we find that temperature is more effective in destroying the spin-glass order than quantum 
fluctuations 

We point out that the renormalisation group trajectory along the critical line flows 
away from the zero temperature fixed point, towards the one controlling the classical finite 
temperature spin-glass transition. This has two main implications: 

(1) The shape of the critical line close to F = is analyticS. In the present case our 
results are consistent with Tc(0) — Tc(F) ~ F^. 

(2) The exponents controlling the transition for finite F are the same as those of the 
classical spin glass transition, except at zero temperature where quantum effects become 
dominant. 

We can develop a scaling theory for the spin glass transition close to the unstable zero 
temperature fixed point at (F/J) ^ (F/J)^. Introducing an exponent z through the scaling 
relationS 

J' = b-'J, (6) 
we obtain the scaling form for the free energy density close to (F/J)^ a 
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X I \2—a ^ 



{T/J) 

m 



(7) 



where g = \ (r/J) — (r/J)c|, i' and z are the correlation length and the dynamic exponents, 
respectively, a is a critical exponent which describes the singularity of the ground state 
energy; all these exponents are associated with the zero temperature fixed point and are 
related through the modified hyperscaling relation 2 — a = v{d+z). Close to (F/ J) ~ (F/ J)^ 
the critical temperature vanishes as 

T, oc \gr (8) 

which allows us to obtain the product vz. Our RG results for the behaviour of the phase 
boundary near T = yield vz ~ 1.23. An independent calculation of the exponent z 
through a finite size scaling analysisil for the gap at the critical point yields z = 1.40. Using 
these results we obtain v = 0.87 consistent with the exact constraint v > 2/d for disordered 
system^. The finite slope of the phase boundary close to T = is a consequence of the fact 
that vz > 1 for the disordered case differently from the pure three-dimensional case where 
vz = 1/2 < 1. 

In conclusion, we have examined the phase diagram of the transverse Ising spin glass 
model. We have compared our data with those obtained experimentally for the dipolar 
Ising spin-glass LiHoo.i67Yo.833F4 with an applied field in the transverse direction. In spite 
of the long-range character of the interactions between spins, this system is more apropri- 
ately described by a model with short range interactions than by the equivalent Sherrington- 
Kirkpatrick model in a transverse fieldi. However the small value of the exponent 7 associ- 
ated with the non linear susceptibility was found experimentally0 to be quite different from 
that of the classical spin glass transition. This remains a puzzle from the point of view of 
our results. Further theoretical and experimental studies are required to clear this point. 
Work is in progress to investigate the possibility of the transition at T = being first order, 
and to obtain a more detailed scaling analysis of this transition. 
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FIGURES 



FIG. 1. Clusters used in the RG transformation in three dimensions. The terminal sites are 
labeled 1 and 6. 



FIG. 2. Critical temperature vs transverse field for the Ising model. (□) and (A) denote 
our renormalisation group (RG) results in the pure and spin glass cases, respectively. The 
replica-symmetry-breaking (RSB, Ref. 9) results (— — — ) and the experimental data (•) for 
LiIIoo.i67Yo.833F4 (Refs. 10 and 11) are also shown. 
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